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Abstract 

We study hysteresis in a micron-sized bead: a non-magnetic matrix embedded with super- 
paramagnetic nanoparticles. These hold tremendous promise in therapeutic applications as heat 
generating machines. The theoretical formulation uses a mean-field theory to account for dipolar 
interactions between the supermoments. The study enables manipulation of heat dissipation by a 
compatible selection of commercially available beads and the frequency / and amplitude h of the 
applied oscillating field in the labortory. We also introduce the possibility of utilizing return point 
memory for gradual heating of a local region. 

PACS numbers: 75.75.Jn,76.60.Es, 75.60.-d 
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The use of small magnetic particles for diagnostic treatments is an emerging area of re- 
search. They provide several exciting opportunities for site-specific drug delivery, improving 
the quality of magnetic resonance imaging, manipulation of cell membranes, etc. UQ. The 
main reason behind this applicability is the ease with which they can be detected and di- 
rected by the application of an external magnetic field. They can be prepared in varying 
sizes, from nanometers to micrometers, comparable to the dimensions of biological enti- 
ties such as genes, proteins, viruses and cells. When coated with relevant compounds, the 
particles can be made to bind to these entities thereby providing a controllable means of 
tagging them. The nanometer-sized particles are usually single-domain and superparamag- 
netic (SPM) in nature. Their relaxation behavior is characterized by Neel and Brownian 
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relaxation times which are strongly size-dependent |4], |5J. The micron-sized particles, com- 
monly called beads, are usually SPM particles packed in a non-magnetic matrix or a porous 
polymer with SPM particles precipitated in the pores . The micron-sized beads have 
several advantages over their nanometer-sized counterparts: they are easy to observe and 
manoeuvre and can be bound to cells due to their compatible dimensions. 

When sufficiently large in number, the supermoments inside the bead start interacting 
via dipole-dipole coupling. They are no longer paramagnetic due to the appearance of a 
dipolar field. The presence of a permanent dipole moment opens up surprising avenues 
in biomedicine. An important non-equilibrium property exhibited by magnetic systems is 
that of hysteresis. Typically, when an oscillating magnetic field is applied, the response 
of the system is delayed leading to hysteresis. The magnitude of hysteresis is determined 
by the competition between experimental time scales (measured by the inverse frequency 
of the applied field) and the relaxation time scale of the magnetic moment. The area of 
the hysteresis loop is a measure of the heat dissipated in the system. The phenomenon 
of hysteresis can hence be advantageously used to heat magnetic beads positioned at a 
given site in the body, perhaps the site of malignancy, to selectively warm it! This effect, 
referred to as hyperthermia in the medical literature, makes it possible to destroy cells or 
introduce a modest rise in temperature to increase the efficiency of chemotherapy. While 
the possibility of using magnetic beads for hyperthermia holds tremendous promise, its 
successful application is a challenge as yet. The reason lies in the several complex questions 
which need to be answered in this context. For instance, what are the clinically acceptable 
values of the field amplitude and frequency of the applied oscillating field to attain maximal 
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heating? Or what should be the heat deposition rate to sustain tissue temperatures of 42° 
C or more? What is the adequate number of particles to be delivered at the specific site for 
achieving the required rise in temperature? It is therefore of great importance to understand 
heat generation and the role of tunable parameters for controlled heating of magnetic beads 
delivered on malignant cells. 

In this letter we first motivate a mean field theory to incorporate the dipolar interactions 
between the magnetic nanoparticles embedded in the bead for the case of large uniaxial 
anisotropy vis-a-vis the thermal energy. The relaxation dynamics can then be described in 
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terms of a two-state rate theory [7j. We study the role of volume fraction of SPM particles 
in the bead and the amplitude and frequency of the applied field on the hysteresis loop area. 
We quantify the dependence of loop area on the above parameters in the form of scaling 
laws. We also discuss the applicability of return point memory for progressive heating of a 
local region. 

In the case of uniaxial anisotropy (in the z-direction say), the magnetic energy of a 
single- domain, SPM particle is given by [?| 



where V is the magnetic volume of a particle, K is the effective magnetic anisotropy constant 
and $ is the angle between the z-axis and direction of the "super" magnetic moment of the 
single-domain particle. Minimum energy occurs at $ = and it and these angles define 
the two equilibrium orientations of the magnetic moment. If VK ^> ksT, where ks is the 
Boltzmann constant and T the absolute temperature, the magnetic moment is mostly locked 
in two minimum energy orientations resulting in the so-called Ising limit. In this limit, 
may be viewed as a dichotomic Markov process in which it jumps at random between the 
angles and it. The jump rate governed by the Arrhenius-Kramer formula 



where A D is the "attempt" frequency. The reciprocal of the jump rate is the Neel relaxation 
time r N . 

In order to incorporate the effects of dipolar interactions, it is essential to add the con- 
tribution due to the "dipolar Hamiltonian" to the magnetic energy of the particle defined 



E = VK sin 2 $, 



(1) 




(2) 
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by Eq. (TO . In the limit of large anisotropy, this contibution is given by 

7-Ld = y^^iljti 2 ^ 1 -j£21^!L ' m zi m zj , (3) 

where 7$ and 7j are the gyromagnetic ratio of i th and j th particle respectively, r~lj is the 
distance between the two particles, is the angle between flj and the anisotropy (z) axis 
and m Z i is the magnetic moment of the ith nanoparticle along the anisotropy axis. Note that 
the angular dependence (1 — 3cos 2 ^j) implies that the interaction can change sign thereby 
switching from ferromagnetic for angles close to the anisotropy axis to antiferromagnetic 
for intermediate angles (55° < % < 125°). Since m is proportional to the volume of the 
particle, Eq. ([3]) can be rewritten as 

H d = jtf^V 2 Ti ^^^^^ os^cos^, (4) 

where /x is the magnetic moment per unit volume and $ is as defined in Eq. (JT]). 

A mean- field approach is invoked to treat the complicated interaction in Eq. ([3]) (9]. 
Each particle is visualized to be experiencing an effective local magnetic field H due to the 
surrounding (magnetic) medium. Further, the Ising limit allows for the replacement of cos$ 
by a two-state variable a. These simplifications yield the mean-field dipolar Hamiltonian: 

nr = nVaH, (5) 

with the field H evaluated self-consistently: 

H = yAV{a) = fxAVtanh I J . (6) 

The parameter A contains all the constants. It fluctuates in sign and has a magnitude 
dependent on the volume fraction of the SPM particles in the bead. 

The presence of H alters Eq. (T2]). Neglecting terms of order H 2 / K 2 due to the assumption 



of large anisotropy, yields the following generalized rate constants 10] : 



V-mt = W- ^±M ), (7) 
A„ = A„exp(- y < A ; B -^> ). (8) 

Kb± 

The relative populations of magnetic moments aligned along and against the anisotropy 
axis (corresponding to $ = and tt) are now described by a master equation for a two-state 



system [9|: 

— n (t) = -A ^ 7r n (t) + X^on^t), (9) 

where hq and denote the fraction of particles with orientation and 7r respectively. 
Solving Eq. (Q, we can obtain the time-dependent magnetization 

m(t) = Vfi[n (t)-n n (t)]. (10) 

We now subject the bead to a time- dependent oscillating field h(t) = h Q cos In ft. It is 
necessary to replace the dipolar mean-field H in the exponent of Eqs. (I7j) and (JHJ) by a time- 
dependent field H(t) = H + h(t). The hysteresis loop can then be obtained by evaluation 
of m(t) using Eq. (ITUl) at each value of the applied field h(t). We allow m(t) to evolve for a 
few cycles of the magnetic field to provide sufficient time for the transients to settle down. 
As the dipolar mean field H is weak, the magnetization reversal is primarily driven by the 
applied field h(t). However we shall see that though weak in strength, it significantly affects 
the relaxation behavior of the magnetic bead. 

It is appropriate to keep a note of the parameter values that are relevent in the context 
of local heating of cells using magnetic nanoparticles. The temperatures required in hyper- 
thermia and chemotherapy treatments are usually in the range of 42°C to 45°C jl|. The 
frequency / of the applied oscillating magnetic field is in the 50 to 1500 KHz range and the 
field amplitude usually varies between 1 to 200 Oe [l|. Several groups have reported that 
exposure to fields with a product h D ■ f not exceeding 4.85xl0 8 A mT 1 s^ 1 is clinically safe 
[l|. We have ensured that our evaluations include the aformentioned ranges of frequency 
and field amplitude. Further, the problem of tissue cooling due to presence of blood flow is 
difficult to address due to its mathematical complexity. However an often-used thumb rule 
is that a heat deposition rate P of 100 mW cm~ 3 is sufficient in most situations 1|. 

Commercially available beads contain nanoparticles of iron, nickel or cobalt compounds 
with tailored volume fractions. Our numerical evaluations of the loop and its area have been 
made for beads comprising of magnetite (Fe^O^) nanoparticles having a diameter of 8 nm 
and an anisotropy constant of 4.68xl0 5 ergs/cm 3 11]. A 1% volume fraction in this case 
corresponds to the presence of approximately 10 4 magnetite nanoparticles in the micron- 
sized bead. For body temperature T = 37° C, the Ising limit VK ^> k^T that we assume 
holds for diameters in the range 6 to 12 nm. We find that our qualitative observations 
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remain unchanged, even in the case of MNP of alternative compounds, as long as we work 
in this limit. 

In Figure 1, we plot the response of the bead to the applied field h(t) for several values of 
the volume fraction Vf. (i) At low volume fractions (Vf ~ 1%), the supermoments continue 
to remain non-interacting and the reversal is governed by the Neel relaxation time. The 
magnetization curve is thus the Langevin function, (ii) As the volume fraction is increased, 
the supermoments begin to interact via dipole-dipole coupling. Our numerics indicate that 
that the resulting mean field H is small for these volume fractions. As a consequence the 
rates Ao^ and A^o governed by Eqs. ([7]) and ([S]) respectively have comparable magnitudes 
and both contribute to the relaxation dynamics, especially near h(t) ~ 0. The approach to 
an all-up or an all-down configuration of the supermoments is thereby impeded, resulting 
in smooth and continuous hysteresis loops, (iii) For high volume fractions (Vf > 30%), the 
mean field H is large enough to make one of the relaxation rates dominate over the other for 
any value of h(t). The system then exhibits a sharp reversal from an all-up to an all-down 
configuration and vice- versa under the action of the oscillating field h(t). 

Figure 2a shows the effect of the field frequency / on the hysteresis loop. Leaf-shaped 
loops are obtained when / ~ A. These go over to elliptic loops with their major axis parallel 
to the h axis on increasing /. For / ^> A, the supermoments are unable to respond to 
the field variation and no hysteresis loops are observed. Figure 2b shows the effect of the 
amplitude of the oscillating field h a on the loops. Their shape ceases to depend on h when 
it is sufficiently large to produce saturated loops. 

In Figure 3, we quantify the effects of Vf, f and h a on the loop area in the form of a 
scaling law. Our numerical evaluations suggest a power law variation A(Vf, f, h a ) ~ V^f^h h 
where the exponents a, f3 and 5 are distinct for cases (ii) and (iii) discussed in the context 
of Figure 1. The values of these exponents are 1.3±0.02, 0.95±0.02 and 1.5±0.02 for case 
(ii) and 2.0±0.02, 0.9±0.02 and 1.1±0.02 for case (iii). The heat deposition rate P can be 
obtained from the loop area A via the relation P — f x A/V where V is the volume of the 
bead. This relationship enables us to identify the clinically useful range of the loop area 
which is indicated in the figure. The scaling law can therefore be utilized to select (several) 
triplets (Vf, f, h Q ) to generate a desired heat deposition rate. Some of them, for the encircled 
values of A which yield P ~ 100 mWcm~ 3 , have been provided in Table I. 

Finally, we study the phenomenon of return point memory in the bead. If the field h(t) is 
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made to cycle with a reduced amplitute to generate a sub-loop, the system returns precisely 
to the same state from which it left the outer loop. This same memory effect extends 
to sub-cycles within cycles. The system thus remembers a heirarchy of states in its past 
external fields. Figure 4 illustrates the return point memory effect in the magnetic bead. 
The inset evaluates the heat dissipated as a function of h Q , the maximal field amplitude 
used to generate the minor loops. We find that return point memory is always observed in 
standard leaf-shaped loops with saturation. This phenomenon, though unexplored in the 
context of hyperthermia, seems a promising technique for gradual heating of a local area. 

In conclusion, the present letter elucidates hysteresis in a micron-sized bead described as 
a non-magnetic matrix embedded with superparamagnetic nanoparticles. The study enables 
manipulation of heat dissipation by a compatible selection of commercially available beads 
characterized by the volume fraction V) of the superparamagnetic nanoparticles and the 
frequency / and amplitude h Q of the applied oscillating field in the laboratory. We hope 
that the simple mean field framework proves useful in efficient design of experiments in 
the context of therapeutic applications such as hyperthermia and chemotherapy. We also 
hope that this study provokes experiments on the utility of return point memory in these 
applications. 
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Figure 1 



Figure 1: Typical hysteresis loops for several values of the volume fraction Vf specified in 
the figure. 
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Figure 2 

Figure 2: Variation in the hysteresis loop as a function of (a) the applied frequency / and 
(b) the amplitude h of the applied oscillating field. 
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Figure 3: A scaling plot which demonstrates that the area of the hysteresis loop scales as 
A(Vf, f, h Q ) oc V^f^h & . The scaling exponents are specied in the text. The region within 
the dotted lines yields heat deposition suitable for hyperthermia. Several triplets (Vf, f, h Q ) 
corresponding to the encircled value of A are specified in Table I. 
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Figure 4 



Figure 4: Hysteresis loop showing return point memory for several minor loops. The inset 
indicates the heat dissipated as a function of h , the maximum amplitude of the applied 
oscillating field for the minor loop. 
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